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Abstract 

A complete suppression of the exponential decay in a qubit (interacting with a squeezed vac- 
uum reservoir) can be achieved by frequent measurements of adequately chosen observables. The 
observables and initial states (Zeno subspace) for which the effect occurs depend on the squeezing 
parameters of the bath. We show these quantum Zeno dynamics to be substantially different for 
selective and non-selective measurements. In either case, the approach to the Zeno limit for a finite 
number of measurements is also studied numerically. The calculation is extended from one to two 
qubits, where we see both Zeno and anti-Zeno effects depending on the initial state. The reason 
for the striking differences with the situation in closed systems is discussed. 
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I. INTRODUCTION 



Measurements modify the dynamics of a quantum system, and the quantum Zeno effect 
(QZE) pQ is the suppression of transitions between quantum states by frequent measurement. 
This result has attracted much attention since it was first discussed. It has been pointed out 
that the opposite effect, i.e., acceleration of a transition, may also occur; this is known as the 
anti-Zeno effect (AZE). A number of experiments claim to have verified both the QZE [2] and 
AZE [3] and some others are planned [3]. The Zeno and anti-Zeno effects are not restricted 
to coherent (reversible) quantum systems: both have been discussed for irreversible decay 
processes [SHH]. Recently, Maniscalco et al. [TU] have investigated the QZE-AZE crossover in 
the model of a damped harmonic oscillator by controlling both the environmental parameters 
and the system-environment coupling. A similar study [IT] was made for a two-level system 
interacting with a spin bath. Somewhat curiously, the most remarkable practical application 
of the QZE consists in reducing (and eventually suppressing) decoherence to control the state 
of a quantum system [6H9]. 

A number of interesting schemes have been proposed in the last few years to counter the 
effects of decoherence. Among these, there are quantum error-correcting codes [12], schemes 
based on feedback or stochastic control [32], quantum optimal control [2], decoherence-free 
subspaces and noiseless subsystems [T3], and mechanisms based on frequent unitary "bang- 
bang" pulses and their generalization, quantum dynamical decoupling [16]. In this context, it 
has been proposed that dynamical decoupling can be unified with the basic ideas underlying 
the QZE [SI |9]. A general discussion of controlling the decay rate using the QZE was given 
in Refs. and [9j here we restrict ourselves to the Zeno dynamics of a two-level system, for 
the case that the measurements probe only the system but not the bath variables. However, 
we do allow for bath squeezing, which is interesting from an experimental viewpoint [TTJ [JB] 
and leads to particularly rich Zeno dynamics: it was shown that measurements sometimes 
reduce the decay but not always [12]. Moreover, by frequent measurements of properly 
chosen observables, a complete suppression (total Zeno effect) of the decay is possible when 
the system is prepared in particular eigenstates of those observables [20]. The choice of 
observables depends on the squeezing parameters of the bath. 

In Refs. [19] and |20l the survival probability was calculated in the limit of high-frequency 
measurement and only selective measurements were considered; we review some results in 
Sec. [IT] However, non-selective measurements have also been found to be relevant in the 
context of quantum coherence control [3 EI]- In Sec. Ill, we show (for the same system) 
that the characteristics of the (anti-)Zeno effect differ significantly between selective and 
non-selective measurement. In Sec. IV, we study the issue in a numerical framework, valid 
for any finite number of measurements. Section [V] contains the extension from the single- 
to the two-qubit case, where we see both Zeno and anti-Zeno effects for multiple states. We 



then discuss our main observations and make some concluding remarks in Sec. VI 



II. ZENO EFFECT: CLOSED AND OPEN SYSTEMS 

Let us consider a closed quantum system evolving under a Hamiltonian H, according to 
the Schrodinger equation 

|w» = m Hm)) ■ (1) 
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We restrict the discussion to a two-level system or qubit, and consider projective measure- 
ment of an observable <t m having the spectral decomposition = \ — IV^KV^ I- The 
index /i represents a particular measurement direction on the qubit Bloch sphere. Suppose 
the system is initialized in eigenstate then the probability of obtaining the eigenvalue 
+1 after evolving for a short time r > is given by 



p(r) 

where 



exp 



iHr 



l-er\ (2) 



e=(W\H 2 \r i )-W\H\r i ) 2 )/h 2 . (3) 

If one considers n successive projective measurements (separated by time steps r) within a 
time interval t = nr, then the probability for the system to remain in the state ) is 

p(nr) ~ (l — er 2 )" . (4) 

For very frequent measurement fi>l, the survival probability is 

p(t) = p(nr) ~ (1 - e(t/n) 2 ) n ~l--t 2 (5) 

which goes to 1 in the limit n — > oo, corresponding to a complete Zeno effect. 

Quantum (anti-)Zeno effects have also been discussed recently for irreversible open sys- 
tems [5] . For a closed system with Hamiltonian dynamics, all states exhibit the Zeno effect 
under frequent projective measurements according to Eq. (J5]), but for an open system, only 
selected states show the QZE depending on the environment. Considering characteristic 
times much longer than the correlation time Tb a th of the bath or reservoir with which it in- 
teracts, the evolution of an open quantum system is well described in terms of the Liouville 
superoperator C by the master equation 

%-C M . (6) 

For a short time interval r, the density operator in terms of its initial value is given by 

p( r ) = e^p(0)^p(0) + £{p(0)}r. (7) 

Taking the system to be a qubit, we next consider measurements of cr^. If the initial state 
is p(0) = IV'iXV'i I; then after time r the survival probability is given by 

V(r) = Tr{p(rM>«|} = WWM) * 1 + «|£{p(0)}|<)r . (8) 

For n consecutive measurements separated by time intervals r, the survival probability for 
the system to remain in the state \ip^) is 

V(nr) ~ [1 + <V£|£{p(0)M>T] n • (9) 

For very frequent measurements and with the total time t = nr kept constant, the 

survival probability becomes 

V(t) = lim V{nr) 



lim 

n— »oo 



l + «|£{p(0)}|<)- 
n 



-i n 



exp[<^|£{p(0)M>*] . (10) 
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It is clear from Eq. (10) that V(t) is in general a function of the initial-state parameters, even 



though this will be left implicit in the notation. Expression (10) is valid only when the time 
r between consecutive measurements is very small, | \£{p(0)}\ijji) | r < 1. At the same 
time, it is assumed that r is nonetheless long enough that (z) the Markov approximation is 
valid for evolution over time r and (ii) the measurement can be taken as instantaneous on 
the time scale r. These together define the regime considered in this paper. 
The freezing of the initial state for frequent measurements is achieved only if 

w\c{ P (o)m) = o . (ii) 

Thus, the Zeno effect in general depends on the dynamics of the system and therefore 
indirectly on the properties of the bath with which it interacts, as well as on the observable 
to be measured. For an open quantum system, it is interesting to ask why only selected 
states show the total Zeno effect unlike the case of closed systems. A key reason is the 
Markov limit [20] already implicit in Eq. (|6]); we shall return to the issue in Sec. VI 



The above will be applied to a two-level system in a broadband squeezed vacuum bath 
(at zero temperature). For this case, the Born-Markov master equation is well studied, and 
given by [EMU M, E3] 

-Qj; = 2^( iV+1 ) ( 2(J -P (J + ~ v+V-P ~ (**+<*-) + ipN (2a + pa_ - (J-O+p - p<7_cr + ) 

- 7 Me i V+pa + - 7 Me- i V_pa_ (12) 
= | (2LpL f - pL ] L - L ] Lp) = £{p} (13) 

in the interaction picture, where L = y/N + 1 o~_ — y/Ne lv cr + , 7 is the decay rate, N = 
sinh 2 (r), M = sinh(r) cosh(r) = a/ iV(iV + 1), with r and 77 being the resonant squeezing 
amplitude and phase of the bath, respectively. 

The effect of this environment is suppressed by the Zeno process for certain states, called 
decoherence-free Zeno states (DFZS). We consider measurement along a unit vector p = 
[sin^) cos(4>), sin(#) sin(</>), cos(#)] on the Bloch sphere, so that the measurement operator 
<7 M is given by 

e» = a.fi = \r 1 )w\-mm- (14) 

The eigenstates of cx^ 1 are 



K) = cos(0/2)|O) + sin(0/2)e^|l) , (15) 
K> = -sin(0/2)|O> +cos(0/2)e^l> , (16) 
|0) and |1) being up and down states respectively [23]. If the system is initialized in the 



state p(0) = IV'i ){^\ \ = \ (1 + c? ■ p), then the Zeno-limit survival probability (10) is 



V(t) = exp[«|£{p(0)}|<>t] = exp[ Cl (0,0)t] • (17) 



Using Eq. (13) for C and (15) for l^i), one obtains 

ci(0, 0) = -|(iV + i) (1 + cos 2 (#)) - | cos(fl) - ^ sin 2 (#) cos(20+p) . (1* 
Note how squeezing M^O breaks the symmetry under 0- rotations [25J. 
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III. NON-SELECTIVE MEASUREMENT 



For dissipative systems, it follows from Eq. (10) that transitions out of the initial state in 



general are possible even in the Zeno limit. Since the system is constantly being monitored, 
such events can be observed if they occur, and one has a choice whether to discard instances 
that contain them. If one does so, the measurement is said to be selective — the situation in 
Sec. |TTJ Retaining all instances regardless of intermediate measurement outcomes is called 
non-selective measurement. For a meaningful comparison, let us also consider the latter 
in detail. In terms of the orthonormal projections Pf = I'^iXV'il 



and P 2 



mm\ 



1 — |, the measurement superoperator P M can be written as 

i 

The evolution of the density matrix after n steps of Zeno dynamics is now given by 

p(t) = p(nr)= (>e^) n p(0). 



(19) 



(20) 



Evaluating Eq. (20) first for n = 1 and again taking p(0) = I^iK^iIj one obtains 

px = p^p(o) = + (i - ^i) mm\ , (21) 

where we used Eq. ^ for the short-time evolution, and with the survival probability = 
1 + cit = 1 — (ipi!, |£{p(0)}|^2 ) r ) with ci as in Eq. (18). Keeping terms up to first order in r, 
the density matrix after the second step becomes 



P2 = p^pi = ^mm\ + (i - w w)m\ , m) 

where ^ 2 = ^ + 9> xCl r + (l-^i)c 2 r, with c 2 = m \C{\if^){if^ Note the physical 
origin of the third term for & 2 . it represents probability which was residing in |?/> 2 ) at time r, 
but which returned to \ipi) during r < t < 2r. This additional (positive) contribution is 
specific to non-selective measurement, which therefore increases the survival probability over 
the selective case: for all t, 

&{t) > V{t) . (23) 

In general, at the (n + l)th step one has & n+ \ = + ^ > n c x T + (1 — ^ n )c 2 T. Hence, for 
very short r (Zeno limit), the survival probability evolves according to 

+ (c 2 - Cl )^ = c 2 , (24) 



dt 



with the solution 



c 2 ( 



7 
2 



9>{t) 



C-2 



Cl 



c 2 - Ci 

,2 



c 2 - Ci 

7 



-(C2-C\)t 



(25) 



(n +0(1 + cos 2 (fl)) - | cos(#) + ^ sin 2 (#) cos(20+r/) . (26) 



In Fig. TJ we plot both survival probabilities: (a) V{9,<p,t) in Eq. (17) and (b) £P(6,(j), t) 
in Eq. (25). In Fig. [T^i, the maxima correspond to ci(9, <p) = [26], leading to a total Zeno 
effect V(t) = 1. For arbitrary iV and r), there are two such maxima [20J: 



cos(6 l 



N-M 



7T 



zly 



iV + M 2(iV + M + l/2) 



•'zl 



(27) 
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FIG. 1. (Color online) Analytic results for the Zeno- limit survival probability for both selective (a) 
and non-selective (b) measurements. Other parameters are iV = 1, 77 = 0, and total time yt = 10. 



and (6*22,022) = (flzi, 0zi+7r). In Fig. [TJd, one observes a total Zeno effect = 1] at two 

locations, and the AZE [&(t) — > 0] for two other points. One can show that c 2 — c\ > for 
any 9 and (ft. Hence, from Eq. (25), the survival probability in the long time limit becomes 



&>{9, 



c 2 



c 2 - Ci 



(28) 



We have the Zeno effect for c\ — and anti-Zeno effect for c 2 = 0. Thus, the Zeno points for 



non-selective measurements are again given by Eq. (27), obtained for selective measurements; 
the agreement is consistent with the inequality (|23[). 

Let us now consider an inverted image (9', (ft') = (tc — 6,7c + 0). One verifies from 
Eqs. Q, Q that c 2 (9',(ft') = -^(9, (ft), so that ^(#',0') = 1 - &{0,<j>) by Eq. @. 
Hence, for non-selective measurements, we see the AZE occurring (Fig. [T]d) at the image 
points (# Z ,0 Z ) = (7r — 9 z ,n + Z ) of the Zeno angles mentioned above [27]. Physically, 
fi(6',(j)') = — £i(9, (ft) so that -H- ) under the inversion, leaving P M invariant. Hence, 
p(0)\e,4> and p(0) 1 6»',0' undergo the same evolution (20); in the stationary limit, p(t) becomes 
independent of p(0) so that = 1 — & expresses Tr{p(t)} = 1. Because of the additional 
symmetry 1— > (ft + 7T, in Fig. [TJo one sees a simple reflection under 9 = tt/2. 



IV. APPROACH TO THE ZENO LIMIT 



The analytic expressions (17), (25) for the survival probabilities V(t) and were 
derived in the Zeno limit, so that we could ignore higher-order terms in r for the evolution ([7]). 
If one also wishes to study the approach to this limit, then our method should be generalized 
to finite r. In this case, the action of e £l ~ is conveniently evaluated by integrating the linear 
system (13) numerically on [0,r]. We will again consider both selective and non-selective 
projective measurements. 
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FIG. 2. (Color online) Selective measurements: numerical survival probability V n for total time 
ryt = 10. The number of measurements n is (a) 1 (b) 5 (c) 30 (d) 1000. The reservoir parameters 
are A^ = 1 and w = 0. 



For the case of selective measurements, one merely has to omit the final linearization step 
in Eq. upon which the finite- r counterpart to Eq. ^ becomes 

V n = (Tr [Pfe^{|^)(^|}]) n = m\e CT m)W\m)T ■ (29) 

In Fig. [2j we plot the n-measurements survival probability V n (9,<fi) for four different n, 
keeping the total evolution time t = nr fixed. One observes a partial Zeno effect for small n; 
as n is increased, V n — > 1 for two points exhibiting a total Zeno effect, each corresponding 
to a DFZS. For other states outside the immediate vicinity of the DFZS, we see a rapid 
decay V n — > indicating AZE behavior. Comparing Figs. [2ji and [T^i, one sees that Piooo 
very closely matches the analytical value (17) for V(t). 

For non-selective measurement, we simply iterate n times the operations of evolution 



and projection in Eq. (20). Clearly, since the projected density matrix involves only one 



(probability) variable as in, say, Eq. (21), a more sophisticated approach is possible but not 
required. Hence, the survival probability is 



Tr 



M>«|} 



(30) 



In Fig. |3j & n (9, 0, t) is again plotted for various n and fixed t. Comparing Figs. [2^, and [3^, 
one sees that trivially V\ = since in this case there are no intermediate measurements. 
For n > 1 however, V n and while still related by V n < & n , differ significantly. In 
particular, for the non-selective case, both the QZE and the AZE are limited to two points 
on the Bloch sphere. This is consistent with the prediction (25), and indeed Figs. [3]i and[l|D 
show the numerical ^1000 and the analytical &(t) to be in excellent agreement. Note also 
how the mirror symmetry with respect to 9 = n/2 emerges for large n: for jt = 10, the 



incoherent Zeno process is in the stationary limit while the dissipative evolution (13) itself 



will reach this limit (recovering the symmetry) only for larger t. We have verified (data not 
shown) that this is indeed the case. 
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FIG. 3. (Color online) Non-selective measurements: numerical survival probability for total 
evolution time jt = 10. The number of measurements n is (a) 1 (b) 5 (c) 10 (d) 1000. The reservoir 
parameters are N = 1 and w = 0. 



V. TWO-QUBIT ZENO DYNAMICS AND ENTANGLEMENT 

Finally, we consider quantum Zeno dynamics for two two-level systems under both selec- 
tive and non-selective measurement. In this context, we can study their effect on quantum 
entanglement, for which both protection [6] and even generation [28] have been shown to be 
possible previously. The quantum master equation for two qubits interacting with a common 
broadband squeezed vacuum reservoir is given by 

dp 
di 



L{p} = -y(N + 1) (2S- P S + - S + S-p - pS+S-) 



+- 7 iV(25 +P S_ 



S-S + p — pS-S + ) 



7Me" ) (2S+pS+ - S+S+p - pS+S + ) 
7Me- ! " (2S_p5_ - S_S_p - pSSJ] 



(31) 



where S± = a\. + o\ are collective raising and lowering operators [29J. Thus, the qubits' 
decoherence dynamics are coupled even in the absence of a direct Hamiltonian interaction 
term between them. Measurement (super)operators generalizing a M of Sec. O and P M of 



Sec. Ill now involve four basis states, chosen as 



= cos(a/2)|00) + sin(a/2)e i/3 |ll) , 
|tf 2 > = -sin(a/2)|00) + cos(a/2)e i/3 |ll) 
|tf 3 ) = cos(5/2)|01) +sin(5/2)e ix |10) , 
|tf 4 ) = -sin(5/2)|01) + cos(5/2)e ix |10) , 



(32) 
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FIG. 4. (Color online) Numerical two-qubit survival probability under very frequent (n = 1000) 
selective measurements, for initial states \^>i(a,/3)) (a) and {^3(6, x)) (b). Other parameters are 
iV = 1, r) = 0, and total time jt = 10. 



where < a, 5 < ir, and -f < /3, X < if- 

The numerical scheme is readily transcribed from the single-qubit case. For selective 
measurements, this involves taking £ \— > L in Eq. (29), with L as in Eq. Q. In Fig.gji, wc 
plot the results for initial state \ipi) i-> |^i); since |\l/2(a,/3)) = |^ r i(a+7r, 0)) gives nothing 
new, in Fig. [4]d we use the antiparallel initial state 1^/3). In the Zeno limit, one can also 
transcribe the analytical result (10). For p(0) = the condition (flTJ becomes 



(^iim^iX*!!}]*!) = Q(a, (3) = -j[2N + 1 + cos(a) - 2Msin(a) cos(/3+^)] = . (33) 



It can be shown that Q(a,/3) < in general and the maximum Q(a,/3) = is attained for 
(3 = —T]. Solving Eq. (33), one obtains a = arccosf— l/(2iV+l)], which can be substituted 
into Eq. J32} to find |*i z ) = (JV|00> + Me- ir >\ll))/y/N 2 + M 2 . If the initial state is |* 3 ), 



then the Zeno condition reads 



<tf3|L{|*3>(tf3|}|*3> = -7(2JV + 1) [sin(5)cos( X ) + 1] = 



(34) 



satisfied for (5 = 7r/2,x 



tt 



which evaluates to |* 3z ) = (|01) - \10))/V2 [6J. These 



solutions of both Eqs. (33) and (34) are in agreement with the respective numerical results 



in Fig. |4j Note that |^i z ) and the maximally entangled singlet |^3 Z ) are decoherence-free 
states of two qubits in a common squeezed bath in the Markovian regime [TBI I22H2I] , and 
thus are also DFZS's in the selective measurement case. 

Having thus verified the two-qubit numerics, let us turn to non-selective measurements. 



In the Zeno dynamics (20) we substitute C t— > L, taking the measurement projector as 
Pp = YliiPiPPi with now Pi = in terms of the states (32). In Fig. [Sa, we plot the 

survival probability ^1000 of the initial state |^/i(a, /?)), with one Zeno and two anti-Zeno 
points visible. In contrast, if the qubits start in \^ 3 (6,x)} (Fig- [5b), one observes a Zeno 
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FIG. 5. (Color online) Numerical two-qubit survival probability under frequent (n = 1000) non- 



selective measurements in the basis (32), for N = 1, rj = 0, and jt = 10. (a) Initial state |vPi(a;, /?)), 
with fixed 6 = n/2 and % — 0- (b) Initial state | ^3 (<5, x)), with fixed a = ir/2 and (3 = 0. 



point but no total AZE. Since Fig. [5] shows only two cross-sections through the parameter 
space for Eq. (32), which in turn presents only a small selection of all possible bases, it is 



thus seen that the two-qubit Zeno dynamics are richer still than in the one-qubit case. 



VI. CONCLUSION 



In summary, we have formulated an accurate numerical scheme to probe the Zeno dynam- 
ics of a qubit evolving under the interaction with a squeezed reservoir. The method is valid 
even for a finite number of intermediate measurements, and coincides with analytic results 
in the limit of infinitely frequent measurement. We have also extended the calculation to 
the two-qubit case, where we see both Zeno and anti-Zeno effects for multiple states. 

For a closed Hamiltonian system, all states show a total Zeno effect as deduced in Eq. 
It is intriguing that for the irreversible system ( 13 ) one finds a total QZE only for selected 



states, since it is supposedly derived from Hamiltonian dynamics for the full system of qubit 
plus bath. Indeed, writing our measurement operator more precisely as S M = <g> Ibatiu 
we are exactly in the situation of Ref. [30j with E M partitioning the total Hilbert space as 
Ti = H\ © %2, where % = spand^f)} <g> "Hbath {i = 1, 2) [21] • Yet, our findings contradict 
Ref. [2D1 even though we do not fault the latter's analysis. Namely, in its Eqs. (l)-(9), Ref. I3"01 
concludes that in the Zeno limit the density matrix components should stay confined to the 
various measurement eigenspaces; for our case, this would imply a total QZE for the qubit, 
with only the bath having some freedom to evolve [32]. The proposed resolution has been 



hinted at below Eq. (11): our starting point (13) is already given in the Born-Markov 



approximation Tbath — > 0, so our "Zeno limit" in fact probes an intermediate regime where r 



is negligible on the scale of qubit decoherence as indicated below Eq. (10), while nonetheless 
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t 3> Tbath- If, however, the time r between consecutive measurements would be so small 
that r <C Tb at h, the dissipative system will exhibit the total QZE, for any measurement 
direction [20]. It would certainly be instructive to pursue the issue further in a unified 
treatment for all r/r bath , but deriving exact non-Mar kovian qubit dynamics in the presence 
of both the squeezed reservoir and frequent measurements is highly nontrivial [33J. 

The cascading time limits above may seem complex, but should be familiar from the 
textbook derivation of Fermi's Golden Rule. The latter describes only the short-time linear 
part of, say, an atom's exponential decay, yet its derivation by time- dependent perturbation 
theory already involved a t — > oo limit (to enforce the resonance condition). 

While the considered regime thus may not be "Zeno dynamics" in the sense of Refs. Q] 
and EIJJ this is not a shortcoming as long as the difference is understood. On the contrary, 
the results (17), (25) and their various generalizations are quite subtle. In particular, only 
because for our model the QZE in general is not total is there such a rich interplay between 
selective and non-selective measurement — a mainly technical distinction for closed systems. 
In view of the proposals for physical realization of squeezed environments [17], it will be 
interesting to see if our predictions can be verified experimentally. 
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